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A self-consistent formalism for the double beta decay of Fermi type is provided. The particle- 
particle channel of the two-body interaction is considered first in the mean field equations and 
then in the QRPA. The resulting approach is called the QRPA with a self-consistent mean field 
(QRPASMF). The mode provided by QRPASMF, does not collapse for any strength of the particle- 
particle interaction. The transition amplitude for double beta decay is almost insensitive to the 
variation of the particle-particle interaction. Comparing it with the result of the standard pnQRPA, 
it is smaller by a factor 6. The prediction for transition amplitude agrees quite well with the exact 
■ result. The present approach is the only one which produces a strong decrease of the amplitude and 

at the same time does not alter the stability of the ground state. 
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I. INTRODUCTION 



< 

CNl , Many body theories have been very successful in describing the collective properties of nuclear systems 11 2[| . Many 
properties are described by considering only the two-body interaction between protons or between neutrons. However 
there are several features which cannot be explained if one ignores the interactions of protons and neutrons [fi3|-p6] . 
. For example the pairing interaction between protons and neutrons seems to be important for determining both the 
ground state and the high spin states properties of nuclei having Z = N [^8| , ^p|j2^j25[] . The monopole and dipole 
proton- neutron interaction was intensively studied in connection with the double beta decay process This 
' subject was considered by many theoreticians especially in connection to the neutrino- less double beta decay (Ov/30). 

Indeed, this phenomenon if it exists, might be crucial for answering a fundamental question of physics, namely whether 
■ the neutrino is a Majorana or a Dirac particle. In order to make predictions in this field one needs reliable tests for the 
nuclear matrix elements. Unfortunately, up to date, such tests are missing. However, since essentially the same matrix 
. elements are used for the two neutrino double beta decay {2vf3(3) for which plenty of data exist, the idea of using 
those NN interactions, which provide a realistic description of the 2vj3j3 process has been widely adopted. The present 
status of the theories proposed for 2vj3j3 is as follows: The approach which produces results closest to experimental 
data is the proton- neutron quasiparticle random phase approximation {pnQRPA). The predictions for the transition 
amplitude is however too high comparing it to the existent data. Therefore one has looked for modifications of the 
existent formalism in order to obtain the desired quenching of the theoretical result. The first new approach was 
proposed by Cha |3l[] for a related subject. Indeed he noticed that the single beta transition rate is very sensitive to 
varying the strength of the two-body interaction in the particle-particle (pp) channel. This idea was adopted by most 
of the groups working in the field of 2vf3(3 p^-|36[|. The results are as follows: The transition amplitude is almost 
insensitive to the variation of the strength, g pp , with which one multiplies the matrix elements of a realistic force, of 
the pp interaction for a large interval starting from zero. Then the amplitude is decreasing quickly reaching zero for 
9pp ~ 1 which, as a matter of fact is very closed to its critical value where the pnQRPA breaks down. It should be 
mentioned that g pp — 1 is just what a realistic force produces. Of course on the fast down sloping part of the plot the 
experimental data is met so one could assume that the problem is solved. However, in this region the pnQRPA is not 
a good approach since the results are not stable against adding new correlations. A positive feature of the pnQRPA 
is that the Ikeda sum rule (ISR) is fully satisfied for any value of g pp . Several attempts have been made to stabilize 
the pnQRPA ground state, or to move the value zero of the transition amplitude to the unrealistic region. Among 
the proposed approaches are the boson expansion technique (BE) |37],|3^] , the multiple commutators method (MCM) 



39|,|40[ , the renormalized pnQRPA 41 and the fully renormalized pnQRPA |42|]. The first three methods succeed to 
shift the zero of the transition amplitude to the unphysical region but the Ikeda sum rule is violated by an amount 
ranging from 15 to 30%. Some complementary features of the boson expansion and re- normalization procedure were 
studied in ref. E3] . Ikeda sum rule is the difference of the total strengths for beta minus and beta plus transitions 
and is influenced by the completeness of the states participating in the process as well as by the correlations involved 
in the ground state. Some groups consider Ikeda sum rule as a signature for the fact that the formalism takes care of 
the Pauli principle [Q. However, even if the Pauli principle is preserved, ISR can be violated if the states involved 
do not have a good particle number. 
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Very recently we started investigating the consistency of the pnQRPA treatment of the pp interaction and the 
structure of the mean field |45| , |46|| . Indeed, without exception all previous formalisms included the particle-particle 
proton-neutron interaction in the pnQRPA approach but ignored this interaction when the mean field was defined. 
In the previous paper we investigated this problem within a consistent approach. We have proved that including 
the particle-particle interaction in the single particle mean field prevents the QRPA to collapse when the interaction 
strength is increased, despite the attractive character of the pp interaction. Moreover the RPA energy is an increasing 
function of the interaction strength. Here we shall call the procedure mentioned above as the QRPA with a self- 
consistent mean field, QRPASMF. 

In the present paper we study the double beta Fermi decay by using the wave functions provided by the QRPASMF 
approach. We are studying the Fermi transition since this has the nice feature that for a single j case, the model is 
exactly solvable and therefore one can test the accuracy of the adopted approximations. The paper is organized as 
follows: In Section 2 we review briefly the results obtained in the previous paper. On this occasion we collect the main 
results needed for the present purposes. In Section 3, the QRPA formalism is written for the quasiparticle operators 
which admit the generalized BCS wave function as vacuum. As we showed in ref. [46] the generalized BCS wave 
function is the static ground state of the system under consideration. The exact treatment is presented in Section 
IV. The results for the four systems involved in the processes considered here, are shown in the tables. The equations 
defining the double beta transition amplitude as well as the ISR are given in Section V. The numerical results are 
discussed in Section VI while the final conclusions are summarized in Section VII. 



II. BRIEF REVIEW OF THE CLASSICAL DESCRIPTION OF THE STATIC GROUND STATE 



Since the present formalism is a natural continuation of our previous work |46|, where the particle- hole (ph) and 



particle-particle (pp) interactions contribute to the structure of the ground state, it is worth summarizing the main 



achievements from ref. 46 1. This will help us to define the context of the present subject as well as to fix the notations 
and conventions. 

We consider a heterogeneous system of nuclcons which move in a spherical shell model mean field and interact 
among themselves in the following manner. Alike nucleons interact through monopole pairing forces while protons 
and neutrons interact by a monopole particle-hole and a monopole particle-particle two-body term. Such a system is 
described by the following many-body Hamiltonian: 

H = ( e r - X r)cl jm C Tjm 

_ ^JL \ "* t t _ ^» \ t t 

j,m;j' .rn' j,m;j' ,m' 

+ 2 x(cpi c 'y)°( c L' c pj')o - 2xi(c^c^)o(c„j'C w v) . (2.1) 

c\^ m {c T j m ) denotes the creation (annihilation) operator of one particle of r(= p,n) type in the spherical shell model 

state \r\nljm) = |tjto). The time reversed state corresponding to \r]nljm) is denoted by |t;jto) = (— ) J ~ m \T;j — 
to). Several authors used this Hamiltonian to describe single and double beta Fermi transitions within a, pnQRPA 
formalism @-||,|§. 

The static ground state of this Hamiltonian is described by the stationary solutions of the time dependent equations 
of motion provided by the variational principle (l^Jl^^] : 

sfwH-i-j^Wdt^O, (2.2) 

with the trial function \*5f) having the following form: 

|*) =$(z p ,z;;z n) z*;z pn ,z; n ) = e T ""e T "e T "|0). (2.3) 
The transformations specified by the operators T are given by: 

jm 
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|0) stands for the particle vacuum state. These transformations depend on the parameters z which are complex 
functions of time. The corresponding complex conjugate functions are denoted by z*. The parameters (z, z*) play the 
role of classical coordinates and conjugate momenta, respectively. The classical coordinates are related to the usual 
coefficients U, V entering the Bogoliubov-Valatin (BV) transformations by the equations: 



Zpj Ppj^- ^ PJ ' i Upj — cos2ppj, 
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(2.7) 



Denoting by U the transformation (2.3) which transforms the bare vacuum into the trial function and by a* the 
images of the particle creation operators through this transformation one easily proves that: 
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The static ground state of the model Hamiltonian corresponds to a minimum value for the classical energy function: 

n = 5> pj - X p )(2j + l)V? Sm + Y^ienj Xn)(2j + 1)V^, 
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(2.9) 



where the following notations have been used: 



^• = (W + ITO 

Vl nj = (0\V nj \ 2 + \Vj\ 2 U 2 pj ), 



off ,nj 
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(2.10) 



The Fermi level energies A p and A„ are determined so that the average number of protons and neutrons are equal to 
Z and N, respectively: 



z = ^(2j + i)y e 2 ffiPj , 

3 

7V = E(2.7 + l)K 2 ff ,„ r 



(2.11) 



As in ref. [46] we restrict our considerations to the single j case. In the paper quoted above we determined the minimum 
points of the classical energy Ti 1 by solving the generalized BCS equations consisting of four equations for the gap 
parameters A p , A„, A pn and /?_ and two constraint equations given by (2.11). The solutions determine the variables 
V p , V n , V which by means of (2.7) and (2.3) provide the static ground state. The next step achieved in ref. [46] was 
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the description of the small oscillations of the system around the static ground state. The procedure is equivalent 
to an RPA formalism. However the semi-classical approach used there requires a certain caution when is applied 
to transition amplitudes. The reason is that while the RPA approximation is a quadratic expansion of the energy 
function around its minimum point, the classical counterpart of transition amplitudes may achieve their minima in 
some points of the phase space, different from the energy minimum. In this case the transition amplitudes are to 
be either expanded around their own minima or around the energy minimum point. In the first case it is difficult 
to perform a consistent quantization of both energies and transition amplitudes while in the second case truncating 
the expansion at its second order might be not sufficient. In order to avoid these complications, here we adopt 
the quasiparticle random phase approximation (QRPA) with respect to the quasiparticle operators having the static 
ground state determined so far, as a vacuum state. The new method will be described in detail in the next Section. 



III. QRPASMF FORMALISM FOR THE SINGLE J CASE 



Inserting the coefficients U and V determined by the pairing equations into Eq. (2.8) one obtains a generalized 
Bogoliubov-Valatin transformation for the quasiparticle operators a[, a 2 . Reversing the BV transformation (2.8) and 
denoting by T the resulting matrix transformation, one obtains: 
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(3.1) 



Making use of this transformation we can easily express the single j restriction of the model Hamiltonian in the 
quasiparticle representation 

H qp = E lia \ai + E 22 a\a 2 + E 12 a\a 2 + E 21 a f 2 a! + AxV^ h V ph - 4xiVj p V pp . (3.2) 

The quasiparticle energies are given analytically in the Appendix A. The two-body terms are written in a factorized 
form. Ignoring the scattering terms, since they do not contribute to the QRPA equations, the factors involved are 
linear combinations of monopole two quasiparticle operators: 



pp 



T ^ T ~i< + ^ T ^2 + ^)4 2 + 2*3*4, 
+T~ 1 T: i A 11 + ^{T~T* nl + T~T* n2 )A 12 + T~T* 2 A 22 , 
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pi nl Fry x pi n 2 p 2 nl' p2 n2 
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The following notations for two quasiparticle operators are used: 

AL = V aL-_aL_, k = 1, 2 

vt 



(3.3) 
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(3.4) 



For a compact writing of the forthcoming equations, it is convenient to introduce the following notations: 

rj = a\ 17 rJ, = a\ 2 , rj = a 22 . (3.5) 

Assuming the quasi-boson approximation for the commutators of the operators r£ and their hermitian conjugate, one 
obtains the following set of linear equations of motion: 



J2 (-AkiT] + Bfcjr,) 



[H,T k 
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(3.6) 
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where the matrices A and B are defined in the Appendix B. The QRPA approximation determines a boson operator 
as a linear combination of the T' and T which describes a harmonic motion of the many-body system: 

k 

[H, C t ] =uC\ (3.8) 
[C,Ct]=l. (3.9) 

Equations (3.8) supply us with the QRPA equations for the phonon amplitudes X and Y: 

A BW ^=.m, (3.io) 



-B -A J \ Y J \Y / 
while the boson condition (3.9) provides the normalization equation: 

]T(|A fc | 2 -|y fe | 2 ) = i. (3.ii) 



IV. EXACT TREATMENT 



Restricting the model Hamiltonian to a single j shell, the resulting many-body Hamiltonian can be exactly treated. 
Indeed, its eigenvalues may be found through a diagonalization procedure in a many-body basis. Since in our 
application we consider a system of 12 neutrons and 4 protons distributed in shells of equal angular momenta (-§■), 
and which may decay through the channels of (3 + ,f3~ and 2v(3(3, the model Hamiltonian is to be considered for the 
following 4 systems with (N,Z) = (13, 3), (12, 4), (11, 5), (10, 6), respectively. In what follows we shall refer to the 
above mentioned systems as to the grand mother, mother, intermediate and daughter nuclei, respectively. To specify 
the paternity of the many-body basis, they will be accompanied by the indices (gm), (m), (i), (d), respectively. Here 
the following basis, formed out of non-orthogonal states, will be used: 

K,n 2 ,n 3 ) =Af nin2n Mi P ) ni (AL) n2 (4n) n3 \0) (4-1) 
where Aj, , A^ n , A^ n are monopole two-particle excitation operators: 



At =Vr t r 1 * 
PP Z_^i Pjm pjm' 
m 

At = V f * J 



4t =Vr t r 



t 

njm 



(4.2) 



The powers ni,n2,ri3 appearing in the defining Eq. (4.1) are different for the four nuclei mentioned above. Indeed 
they are subject to the restrictions: 

2ni + n 3 = Z, 2n 2 +n 3 = N. (4.3) 

The solutions of the above equations for integer numbers define the basis states for the nuclei involved in the transitions 
considered in the present paper. These are given in Table 1. 

The states are normalized to unity. The normalization factors, denoted by Af ni n 2 n 3 , are listed in Table 2 while the 
overlaps of states are collected in Table 3. 

In order to find the eigenvalues of the model Hamiltonian H in the non-orthogonal basis we have to calculate its 
matrix elements. This goal is achieved in Appendix C. Let us denote by H the resulting matrix associated to H. Also 
we denote by O the overlap matrix. As we showed in ref. [46], the eigenvalue equation: 

H\$)i = / = gm, m, i, d, (4.4) 

with 
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|d>), = £c fe |fc>* (4.5) 

k 

can be reduced to an ordinary eigenvalue problem for a symmetric matrix: 

HX = EX, (4.6) 

where 

H = U- 1 H{U- 1 ) T , X = U T C, (4.7) 
while U is determined by the factorization: 

O = UU T . (4.8) 

Therefore the eigenvalues E arc determined by a standard diagonalization procedure for the symmetric matrix H and 
the eigenfunction |$) is determined by the column vector 

C=(U T )~ 1 X. (4.9) 

These results will be used in the next Section to calculate the transition amplitude for the double beta Fermi transition. 
The dependence of the energies obtained for the mother, intermediate and daughter nuclei, on the strength parameter 
of the pp interaction is shown in Fig. 1. 

V. DOUBLE BETA TRANSITION. 

The double beta decay with two neutrinos in the final state is considered to consist of two consecutive single f}~ 
decays. The intermediate state, reached by the first /3~ transition, consists of an odd-odd nucleus in a pn excited 
state, one electron and one anti-neutrino. If in the intermediate state, the total lepton energy is approximated by 
the sum of the electron rest mass and half of the Q- value of the double beta decay process, the inverse of the process 
half- life can be factorized as follows: 

(Tfyj- 1 = F\M F \\ (5.1) 

where F is a lepton phase integral while the second factor is determined by the states characterizing the nuclei involved 
in the process and is given by the expression: 

M v m (o + \p + \o + k) mm (oi\oi,} dd (ot,\p + \o + ) d 

F tt> E k + AE ( - ' 

where the transition operator is 

+ = E c U<**» (5-3) 

m 

The initial and final states are the ground states of the mother (|0 1 ') m ) and daughter (|0^)d) nuclei. In the QRPASMF 
formalism these are described by the phonon operator vacuum state while in the exact treatment these are the lowest 
states obtained through diagonalization, respectively. The intermediate states are one phonon states in the QRPASMF 
approach while in the exact treatment they are eigenstates of the model Hamiltonian describing the odd-odd system. 
The corresponding eigenvalues, normalized to the ground state energy of the mother nucleus, are denoted by E k . Since 
the excited states corresponding to the two vacua are not orthogonal to each other the overlap matrix elements should 
be inserted between the matrix elements describing the two legs of the double beta transition. For the exact treatment 
the overlap of the eigenstates is obtained by multiplying the vectors of the weights produced by diagonalization with 
the overlap matrix. In the QRPASMF treatment, the overlap matrix elements are: 

/n+ |n+ \ V m yd , ym yd , ym yd ymyd ymyd ymyd (K A\ 

The matrix elements describing single f3 + and /3~ transitions of the mother system satisfy the Ikeda sum rule: 

(3 ( - } - /3 (+) = N - Z, (5.5) 

where and denotes the total strength for f3~ and (3 + transitions respectively, while N, Z specifies the number 
of neutrons and protons. 
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VI. NUMERICAL RESULTS 



We applied the formalism described in the previous sections to a system of 12 neutrons and 4 protons moving in a 
shell with j ' = We study the decay of this system to the daughter nucleus with (TV, Z) = (10, 6). This transition 
is achieved in two steps, with the intermediate odd-odd system having (N,Z) = (11,5). The static solution for the 
equations of motion were found using the following set of parameters for the model Hamiltonian restricted to the case 
of a single shell: 

e p = e n = ZMeV, G p = 0.25MeV, G n = O.UMeV, X = 0-20MeV. (6.1) 

The strength of the pp interaction, xi, is varied freely in the interval [0.0,1.5]MeV. For each value of xi, we solved 
the generalized BCS equations for both mother and daughter nuclei. As a result one obtains the gaps and chemical 
potentials and subsequently the U and V coefficients. In the next step, the QRPA equation is solved. One should 
remark that, as shown in ref. [46], there are two constants of motion which results in having two spurious solutions. 
Therefore there exists only one physical solution which corresponds to the matrix elements An and B n , i.e. the 
degree of freedom accounted by the amplitudes Xi,Yi. The RPA energies and phonon amplitudes are used then to 
calculate the amplitude for the double beta decay, given by eq.(5.2). The results for Mp are given, as function of 
Xi, in Fig. 2. The present QRPA approach is different from the standard pnQRPA, where the phonon operator 
is a neutron-hole proton-particle excitation. Here, due to the proton-neutron mixing in the quasiparticle operator 
O-Xjmi the phonon operator comprises also pn scattering terms (ejejj as well as charge conserving excitation operators 
(c p Cpi). Therefore one expects to obtain a description which is substantially different from the one provided by 
the standard pnQRPA. The double beta transition amplitude has been alternatively calculated by using the exact 
eigenstates of the model Hamiltonian H. Thus, H (2.1) was diagonalized for the three systems involved in the process, 
i.e. the mother, intermediate and daughter nuclei. The corresponding energies are plotted in Fig. 1 as function of 
Xi, for vanishing chemical potentials. Several remarks are worth to be mentioned. Note that due to the specific 
model space, j ' = -y, and the proton and neutron numbers the mother and intermediate nuclei, have three different 
+ states while the daughter nucleus exhibits four independent states. The ground state of the mother nucleus gets 
higher in energy than the ground state of the intermediate nucleus. This happens for xi equal to about 0.9 MeV. 
It is not clear if the zero of the double beta transition amplitude in the standard pnQRPA approach, at this value 
of xi, is caused to some extent by this level crossing. The single beta transition to the odd-odd system is therefore 
virtual up to the intersection point and real from there on. As shown in Fig. 2, for the value of xi where the ground 
states of the two neighboring nuclei, mother and intermediate, are equal to each other, the transition amplitude has 
a minimum value. Increasing further the strength xi the amplitude is slowly increasing and finally reaches a plateau 
at about xi = 1-6. The agreement between the QRPASMF result and the exact one is quite good. The amplitude 
predicted by the QRPASMF formalism is varying very slowly with xi which, in fact, reflects the stability of the mean 
field with respect to the particle-particle interaction. Note that within QRPASMF the Fermi sea energies X p and A n 
depend on xij due to the adopted variational procedure. The exact result shown in Fig. 2, corresponds to those 
A P ,A„ provided by QRPASMF for xi = 1- We checked that the Fermi energies are only slightly depending on xi 
and such a dependence does not essentially affect the exact results. 

The previous theoretical investigations mostly focused on the decreasing part of the double beta decay amplitude 
within the standard pnQRPA description, since the stable result is usually larger than the experimental one. Unfor- 
tunately the predictions in this region are not very stable and moreover show a zero for a certain xi- In this context 
it is remarkable that the present formalism produces the desired moderate suppression of the transition amplitude 
and moreover the result is stable against adding more correlations. 

The T=l proton-neutron interaction was also considered in refs. [23] and [49] to define a generalized pairing field. 
The formalism was employed for the description of the Gamow- Teller 2vpj3 process. However the particle-particle 
two body interaction, whose strength causes the ground state instability is of dipole-dipole type and not included in 
the mean field. Accounting also for this dipole-dipole proton-neutron interaction requires a simultaneous treatment 
of the T=l and T=0 proton-neutron pairing interactions. Note that such a problem does not appear for the Fermi 
double beta transition. Indeed here there is no particle- particle interaction term which is left out when the mean 
field is defined. Extension to the more complex situation of the Gamow- Teller double beta decay is under work and 
the results will be published in a forthcoming paper. 

Let us compare the present results with those given by the standard pnQRPA: We recall that within the standard 
pnQRPA, used widely in the literature, first one treats the proton-proton and neutron-neutron pairing interactions by 
the usual BCS procedure (apart from the work of the Tuebingen group, which, as we mentioned before, included also 
proton- neutron pairing j23|,^9| . In the quasiparticle representation, after ignoring the scattering terms the Hamiltonian 
reads: 

H = E P N. P + E n N n + XtA+A + \ 2 {A+ 2 + A 2 ), (6.2) 
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where N p and N n denote the proton and neutron quasiparticle number operators, respectively. The quasiparticle 
energies E p and E n have simple expressions: 

g n 

E T = — T — , r=p,n. (6.3) 
The proton neutron two quasiparticle monopole operator is denoted by: 

At = -^Y a L m at (6-4) 



/20 

The coefficients Ai and A2 have simple expressions in terms of particle-hole and particle-particle interaction strengths. 

X 1 = 2 x {U 2 p V 2 + UX 2 ) - 2 Xi (U 2 p U 2 n + V 2 V 2 ), (6.5) 
X 2 = 2{ X + Xi)U P U n V p V n . 

The pnQRPA energy can be expressed in terms of the strengths parameters Ai, A 2 , as 



+ E n + AO 2 - 4X1 (6-6) 



The transition amplitude (5.2) corresponding to this mode is represented in Fig. 3 as a function of xi an d will be 
hereafter referred as the standard pnQRPA transition amplitude. Comparing the curves for Mp given in Fig. 2 and 
Fig. 3, one notices that the present formalism predicts an amplitude which is smaller by a factor 6 than the amplitude 
at the stable plateau at small xij given by the standard pnQRPA. Since in the three cases considered here, the 
exact treatment, the new QRPA and the standard pnQRPA, we consider the same energy shift (AE = 2.5MeV) 
in the denominator of eq. (5.2) we have to pay attention to possible different zero point energies in the mother and 
daughter nuclei. What is the reason for such a dramatic change in the magnitude of the transition amplitude? To 
see that, we plotted in Fig. 4 the ratios of the matrix elements describing the two legs of the double beta transition 
(see eq. (5.2)), obtained by the standard pnQRPA and the QRPASMF approach. From there one sees that while the 
0~ decay matrix of the first leg elements predicted by the two approaches are about the same, the matrix elements 
describing the second decay differ from each other by a factor of about six, in the beginning of interval. Although the 
matrix element describing the second decay in the present formalism is small it is only slightly changed by increasing 
the parameter xi- Since the static ground state achieves the minimum for the classical energy one expects that the 
present RPA approach is a very good many-body approach. This is confirmed by the results shown in Fig. 5 where 
indeed the phonon backward amplitude Y does not exceed 10% of the phonon forward amplitude in a large interval of 
Xi- By contradistinction within the standard pnQRPA approach the Y amplitude for the daughter nucleus becomes 
equal to the X amplitude (see Fig. 6), but of opposite sign, at Xi about 1. For the mother nucleus the equality for 
the magnitudes of the two amplitudes is reached at xi = 1-2- At the points mentioned above the pnQRPA phonons 
collapse, respectively. The relationship of the two amplitudes suggests that the phonon operators become hcrmitian 
for this strength of the pp interaction, which prevents a well defined boson excitation. Note that this bad behavior 
does not show up in our formalism. The reason is that the pp interaction is considered in the definition of the mean 
field. The Ikeda sum rule ( ISR) is shown in Fig. 7. We notice that the predicted ISR differs from the N — Z value 
by 25% for xi = 1- Of course the exact result agrees perfectly with the ISR. A possible reason for the discrepancy 
quoted above could be that in the present approach, the mother nucleus can decay by 0~ and by + only to one state 
in each case. By contrast in the exact treatment the ground state of the mother nucleus is linked by non-vanishing 
matrix elements of 0~ to three states in the intermediate odd-odd nucleus and by + to two states in the so called 
grandmother nucleus. 

Before closing this Section we would like to comment the comparison of the three methods, exact, QRPASMF, 
pnQRPA, from a different view angle. Note that the QRPASMF result is depending very little on xi an d, in a way, is 
averaging the behaviour of the exact result around the value of xi where the relative position of the ground states of 
mother and intermediate nuclei is changed. The monotonic behaviour of Mp is consistent with the QRPASMF energy 
represented as function of xi- The common feature, monotony, is caused by the fact that for each value of xi °ne 
determines the corresponding static ground state. Despite the fact the particle-particle proton-neutron interaction 
is attractive the QRPASMF mode does not collapse since, at the same time, the interaction increases the pairing 
energy gap and in this way the mode energy becomes an increasing function. In order to get a deeper understanding 
of the mechanism causing the fact that the Mp value predicted by QRPASMF does not exhibit the allure of the 
similar function produced by the exact calculation, some additional comments are necessary. Indeed, from Fig. 1 it 
results that for a certain value of Xi the ground states for the mother and intermediate nuclei have the same energy. 
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However we cannot speak about a ground state degeneracy which might be related to a phase transition for the 
nuclear system since we deal with two distinct systems characterized by different number of protons and neutrons in 
the framework of the exact treatment. Note that this level crossing point produces a minimum value for Mp. On the 
other hand in the standard pnQRPA formalism, where one treats first the proton-proton and neutron-neutron pairing 
interactions, the small oscillations are treated with respect to a static ground state which is independent of xi- Given 
the attractive character of the particle-particle interaction the mode energy is continuously decreasing and fatally 
reaches the value zero. At this point the ground state of the daughter nucleus becomes degenerate. This reflects 
the fact that the static ground state has two dominant components of the same magnitude. In contradistinction 
to what happens in the exact calculations, the components whose amplitudes become of equal magnitudes are both 
associated to even-even systems. Moreover, in the non-exact treatment the two degenerate states characterise the 
same superfluid system. In this respect the monotonic decreasing behaviours of the M F functions corresponding to the 
exact and the pnQRPA descriptions are caused by different physical circumstances. In the present formalism, when 
one passes from the mother to the intermediate nucleus by subtracting a neutron and adding a proton one obtains a 
state of two quasiparticles. Due to the presence of an energy gap caused by the proton-neutron pairing interaction 
the two systems, the mother and intermediate odd-odd nuclei, will never have equal energies. Moreover the above 
mentioned gap is an increasing function of \i which implies an overall repulsive character of the xi-interaction with 
respect to this type of excitation. Aiming at removing possible confusions we stress the fact that the "exact solution" 
in the present paper has different meaning than in all previous publications (see for example Refs [44,48]). Indeed 
our exact solutions are eigenstates of the many body Hamiltonian (2.1) while all previous publications deal with the 
exact eigenstates of the quasiparticle Hamiltonian (6.2) which is a severely truncated image of the initial Hamiltonian 
(2.1) through the Bogoliubov Valatin transformation. 

Concerning the extension of the formalism to a realistic interaction and a large model space for single particle states 
the amount of work involved can be evaluated as follows. For the Fermi transition the extension is straightforward and 
no principle difficulties appear. As for the Gamow- Teller transition there are two ways to cope with this problem, a) 
To treat simultaneously the T=l and T=0 pairing and project out the good particle number, angular momentum and 
isospin from the resulting generalized BCS wave function. The projected states are to be used in a time dependent 
variational equation to define the QRPA states, b) A deformed single particle mean field can be defined, as usual, 
by the spherical shell model term and the particle-particle proton-neutron two body interaction. At the second step 
the pairing interaction for the deformed single particle states is treated. Finally the usual QRPA procedure accounts 
for quasiparticle correlations. We would like to mention the fact that there is work in progress on these lines. The 
two extensions mentioned above have in common with the formalism of the present paper, the idea that first of all an 
optimal mean field, which assures a stable ground state for the QRPA vibrations, should be constructed by involving 
all types of particle-particle like interaction. 

VII. CONCLUSIONS 

In the previous sections we developed a new formalism for the double beta decay. The QRPASMF formalism formu- 
lated here is different from the standard pnQRPA approach in that the static ground state includes the correlations 
coming from the ph and pp two-body interaction. In the previous work we showed that this defines a stable ground 
state with respect to the RPA excitation. Here we show that the new formalism provides a double beta transition am- 
plitude which does not vanish in a large interval of the pp interaction strength. Moreover this amplitude is diminished 
with respect to that predicted by the pnQRPA in the stability interval of xi, by a factor of 6. Such a suppression 
is due to the fact that the second leg of the double beta transition is very much quenched otherwise the first leg 
being only slightly modified. As a matter of fact this is an important result since it improves the agreement with the 
experimental data for the half life by a factor of 36. Therefore the present formalism produces a moderate suppres- 
sion of the double beta transition amplitude. Moreover the predicted amplitude is not very insensitive to increasing 
the strength xi of the pp interaction. This seems at variance with previous theoretical studies which considered 
the particle-particle interaction very important for describing quantitatively the decay rate in pnQRPA approach, 
because it was not included in the self-consistent mean field. The predicted result for the transition amplitude agrees 
quite well to the exact result. The exact transition amplitude exhibits a minimum, but does not vanish, for the value 
of the parameter xi where the ground state energy of the odd-odd system gets lower than the ground state of the 
mother nucleus. It is an open question whether this feature persists also in realistic calculations and contributes to 
the vanishing of the transition amplitude in the pnQRPA formalism. 

The Ikeda sum rule deviates from the N-Z value by an amount less than 4% for xi — and 25% at xi — 1- A 
possible explanation for this discrepancy is that in the single j level, the QRPA ground state of the mother nucleus 
is (3~ decaying to only one state in the odd-odd nucleus. This feature contrasts to the exact calculations where there 
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are three channels for the beta minus and two channels for the beta plus decay. 

Concluding the results of the present paper constitute a very important test of our idea that the suppression of 
the transition amplitude with increasing strength \i °f the particle-particle interaction is caused by a instability 
of the ground state. This un-wanted feature can be cured by considering the contribution of the particle-hole and 
particle-particle two-body interaction already in the single particle mean field. The agreement of the present approach 
with the exact results are encouraging us to apply the proposed formalism to a realistic two-body interaction and an 
extended model space for the single particle motion. 

VIII. APPENDIX A 

Here we give the analytical expressions for the quasiparticle energies appearing in Eq. (3.2). 

£11 = (e P - X P )(\T pl \ 2 - |T~| 2 ) + (e„ - A„)(|T nl | 2 - |T~| 2 ) + A p T pl T~ + A n T nl T~ (8.1) 
7=(T P iT„i - T~T~) + —^{T vl T~ + T~T„i), 



E 22 = (e P - Ap)(|T p2 | 2 - |T~| 2 ) + (e„ - A„)(|T n2 | 2 - \T~ | 2 ) + A p T p2 T~ + A n T n2 T~ 

E12 = (e p - X p )(T pl Tp 2 - T~T~) + (e„ - \ n ){T nl T n2 - T~{F~ 2 ) 
+ \\(T pl T~ + T p2 T~) + ^A„(T nl T~ 2 + T n2 T~) 

+ ^^ iyTplTn2 + Tp2Tnl ~ T pl T n2 ~ T P 2 T nl) 

E 2 i = E\ 2 . 



IX. APPENDIX B 

By elementary calculations one finds the explicit expression for the matrix elements An- and Bik involved in the 
linearized equations of motion (3.6). They are listed below: 

A n - 2£u + A X (T 2 ~T 2 lX + T 2 pl T\) - ±Xi(T 2 ~T 2 ~ + 1*1*), 

A 12 = V2E 21 + 2V2 X (T 2 1 T~T~ + T^T-T- + T 2 ~T nl T n2 + T 2 ~T pl T p2 ) 

-2V2 Xl (T 2 ~T~T~ + TlT~T~ + T 2 pl T nl T n2 + T^T pl T p2 ), 
A 13 = 4 X (T~T~T nl T n2 + T pl T p2 T~T~) 

- 4 X i(T pl T nl T p2 T n2 + T-T-T-T-), 
A 21 = A 12 , 

A 22 = E 11 + E 22 + 2 X [(T~T nl + T~T n2 ) 2 + (T~T pl + T~T p2 ) 2 } 

- 2 Xi[(T~ 2 T~ + T~T~ f + (T n2 T pl + T nl T p2 )% 
A 23 = V2E 21 + 2V2 X {T 2 ~T nl T n2 + T 2 ~T pl T p2 + T 2 2 T~T~ + T 2 2 T~T~) 

-2V2 X1 (T 2 ~T~T~ + T 2 ~T~T~ + T 2 2 T pl T p2 + T 2 2 T nl T n2 ), 
A 31 = A 13 , 
A 32 = A 23 , 

A 33 = 2E 22 + A X {T 2 ~T 2 2 + T 2 p2 T 2 ~) - A Xl {T 2 ~T 2 ~ + T 2 p2 T 2 2 ). (9.1) 
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Bu - 8(x + X i)T~T~T pl T nl , 

B 12 = 2V2( X + Xl )(T~T~T pl T nl + T~T~T p2 T nl + T~T~T pl T nl + T~T~T pl T n2 ), 
Bis = A X {T~T~T p2 T nl + T~T~T pl T n2 ) + 4 Xl (T~T~T n2 T p2 + T~T~T nl T pl ), 
B 2 \ = B\ 2 

B 22 = 2 X ((T~T nl + T~T n2 )(T~T pl + T~T p2 ) + (T~T pl + T~T p2 )(T~T nl + T~T n2 )) 

+ 2 Xi ((T~T~ + T~T~)(T n2 T pl + T nl T p2 ) + (T n2 T pl + T nl T p2 )(T~T~ + T~T~)) , 

B 23 = 2V2( X + Xi) 

( T i2 T n2( T ^ P 2 + T pl T n2 ) + T p2 T n2 (T~T~ + T~T~) , 

B31 = -B13, 

-S32 = -B23, 

£33 - 8(x + X i)T~T~ 2 T p2 T n2 . (9.2) 



X. APPENDIX C 



To calculate the matrix elements of H in the basis (4.1) we need to know the result of acting with H on a given 
representative state. To give an example let us consider the proton-proton pairing interaction and determine the 
vector expansion: 

Q 

~ ~£ A ip A Pp\ l ) k = ^2'Pim\m} kl k = gm,m,i,d. (10.1) 

m 

Once the expansion coefficients are known, the matrix elements of the pairing interaction are readily obtained: 

- ^{ n \A\ p A pp \l) k = 0^lV lm . (10.2) 

m 

The expansion coefficients associated to various two-body terms of the model Hamiltonian considered for a single 
j-shell are listed in tables 4-7, for the grand mother, mother, intermediate and daughter systems respectively. Using 
the expansion coefficients and the overlap matrix we derived analytical expressions for the matrices H corresponding 
to the basis labeled by gm, m, i, d, respectively. Further, these matrices are treated as described in Section IV. 
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TABLE I. The basis states describing the nuclei "grand mother" \gm), mother \m), intermediate \i) and daughter |d) in a 



single j shell. 





(N,Z)=(13,3) 
i=gm 


(N,Z)=(12,4) 
x—m 


(N,Z)=(11,5) 
x=i 


(N,Z)=(10,6) 
x=d 


Il>- 


10,5,3} 


|2,6,0) 


|2,5,1) 


|3, 5,0) 


\2) x 


|1,6,1) 


|1,5,2) 


|1,4,3) 


|2,4,2) 


|3>* 




10,4,4} 


10,3,5} 


|1,3,4) 


|4>, 








|0,2,6) 



TABLE II. The values of (Af^)' 1 /a with a = (30 x 10!)'. The norm of the state \k) x is denoted by Afjf\ 

labels the + states in grand mother (gm), mother (m), intermediate (i) and daughter (d) nuclei, respectively 

x \l) x \2) x ^ \3) x |4) a 

gm 2 7 x 3 2 8 x 3 x V2 

m 2 9 x 3 x y/Z 2 7 x v^3 

i 2 9 x V3 2 7 x ,/f 

d 2 10 x 3 3i x ^58 



|3>- 

2 6 x 3 x 
2 4 x 

x vTl22 



TABLE III. The overlaps of states characterizing the grand mother (gm), mother (m), intermediate odd-odd (i) and daughter 
(d) nuclei, 
[h] 



Nucleus 



|1> 
1 

s/2 



|2> 

3 
1 



1 3) 



— i — 

3v Tl 

°V 143 



|4> 



1 

3_ 

1 

3vTT 




1 

5 

■v/858 



1 /_3_ 

6 V 143 




3 





34 
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Bp n B pn \l)gm 36 
-'pn -* 

-Ap n Apn | l)gm 24 



^-pn ^-pn 1 2) gm 24 ^. (gm) 



TTTTTTT 



TABLE IV. The actions of various terms of the model Hamiltonian on the states describing the grand mother nucleus 

(N,Z)=(13,3). 

Q\n) gm \l) gm |2) 6 

Ap P App\l} gm —6 , (m 
t 

AppApp\'2}g m 36 

A nn A„„\V) gm 60 — 6 .(gm) 

t 

-^titi 1 2} m 96 



B pn B pn \2} gm 24 ^ gm) 18 
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TABLE V. The actions of various terms of the model Hamiltonian on the states describing the mother nucleus (N,Z)=(12,4). 



Q\n) m 

A nn A nn | l)r 
A nn A nn 1 2} r 

A nn A nn 1 3} r 

^pn Bp n 1 1)77 

^pn ^pn 1 2} 77 

Bp n Bp n (3)77 

Ap n Ap n J l)rr 

Ap n Apn 1 2} 77 
^pn^pi 1 3} rr 



|l)m 

72 

JV. 



(m) 
2 



120 



(m) 
2 



(m) 
2 



|2>, 



|3>, 



32 



-12 



A/V 



60 



-12 



A/V 



18 



A/-J 

30 

(m) 



-48 



14 



36 



20 



2 



40 



-20 



AC' 

36 



— 



TABLE VI. The actions of various terms of the model Hamiltonian on the states describing the intermediate odd-odd nucleus 
(N,Z)=(11,5). 



Q\n)j 



A^ppApp\ l)i 

-^pp^ppl^}-! 

A nn A nn J l)i 

A nri I 2} j 
AnnA nn 1 3} % 
-Bp n Bp n I l)i 
Bp n Bp n |2)j 
^pn ^pn 
^prt^pJi I l)i 
Ap n Apn |2}i 



12), 



13), 



64 



100 

Af (l) 



20 

Af (i) 



28 



-20 



A' 



(0 



AA<'» 



-20 



64 

AT<'> 



40 



38 

Af W 

24 



36 



16 



Af> 
A? 1 

40 



-16 



JV„« 



A/; 

50 
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TABLE VII. The actions of various terms of the model Hamiltonian on the states describing the daughter nucleus 
(N,Z)=(10,6) 



0\n) a 



Ap p App\l) d 

AppAp P \2) d 
AppAp P \S}d 
AppA P p\4) d 

An n A nn 1 1) d 
-A-lmAnn \2) d 

A'l in A rin \3)d 

An n A nn 
Bpn B pn \l)d 
-Bpn -Bpn |2}d 
Bpn Bpn 1 3} d 

B Pn B pn \4) d 

Ap n Apn | l)d 
Apn Apn 1 2} d 
Ap n Apn\3) d 

A\ n A pn \A) d 



|2} d 



|3> d 



14). 



96 



-2 



120 



Af {a > 



-12 



56 



Af< d > 

-12-^br 

Af< d) 



60 



12 



32 



at; 



-60 



AT, 



(d) 



14 



-30 



24 

AT| d » 



AC' 



48 



-30 



32 



/v 



42 



30 



A/i 



(d) 



-32 



(d) 



36 



24 



12 



(d) 



36 



-12 



a-; 



(d) 



Af< d > 



66 
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FIG. 1. The exact energies for the states of the mother (dashed line) intermediate (dashed dot line) and daughter (full line) 
are plotted as function of the particle-particle interaction strength, \i- 
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FIG. 2. The transition amplitudes for the double beta (2v/3/3) Fermi transition, corresponding to the exact eigenstates and 
energies of the model Hamiltonian (full line) and the new QRPA approach of the present paper(dashed line), are plotted as 
function of \i- 
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FIG. 3. The transition amplitude for the double beta Fermi transition calculated within the standard pnQRPA, is represented 
as function of \i- 
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FIG. 4. The ratio of the first single beta matrix element involved in eq (5.2) obtained within the standard pnQRPA and 
the present approach respectively (full line) , is represented as function of xi • The ratio corresponding to the second leg of the 
double beta transition is plotted by a dashed line. 
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FIG. 5. The magnitude of the ratio of the phonon amplitudes calculated within the present approach, is represented as 
function of \i for mother (full line) and daughter nuclei (dashed line). 
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FIG. 6. The magnitude of the ratio of the phonon amplitudes, calculated within the standard pnQRPA, is represented as 
function of xi for mother (full line) and daughter nuclei (dashed line). 
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FIG. 7. The N — Z sum rule corresponding to the exact description of nuclear states(dashed line) an to the present approach, 
is represented as function of \i- 
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